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Abstract 



Exact solutions of Dirac equation in two spatial dimensions in the Coulomb 
field are obtained. Equation which determines the so-called critical charge of 
the Coulomb field is derived and solved for a simple model. 



I. INTRODUCTION 



Planar nonrelativistic electron systems in a uniform magnetic field are fundamental quan- 
tum systems which have provided insights into many novel phenomena, such as the quantum 
Hall effect and the theory of anyons, particles obeying fractional statistics . On the other 
hand, planar electron systems with energy spectrum described by the Dirac Hamiltonian 
have also been studied as field-theoretic models for the quantum Hall effect and anyon the- 
ory |3|H] . Related to these field-theoretic models are the recent interesting studies regarding 
the instability of the naive vacuum and spontaneous magnetization in (2+1) -dimensional 
quantum electrodynamics (QED), which is induced by a bare Chern-Simons term ||. In 
view of these developments, it is essential to have a better understanding of the properties 
of planar Dirac particles in the presence of external electromagnetic fields. 

In this paper we would like to consider solutions of Dirac equation in two spatial di- 
mensions in the presence of a strong Coulomb field, and to discuss instability of the Dirac 
vacuum in a regulated strong Coulomb field. In three space dimensions the effect of positron 
production by strong Coulomb field was predicted in M and were studied in pH[5]. 



II. MOTION OF AN ELECTRON IN THE COULOMB FIELD 

Let us consider a relativistic electron in two spatial dimensions in a Coulomb field the 
vector potential of which is specified as 

A°(r) = -Ze/r, A x = A y = 0. (1) 

In 2+1 dimensions the Dirac matrices may be represented in terms of the Pauli matrices. 
We choose the representation a = (—a 2 , a 1 ) and (3 = a 3 . Then the Dirac equation has the 
form (c — h—1) 

(idt-H D )* = 0, (2) 

where 
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H D = a ■ P + 0m + eA° = c^Pa - a 2 Pi + cx 3 m + eA c 



(3) 



is the Dirac Hamiltonian, = id^ — eA^ is the operator of generalized momentum of 
electron, m is the rest mass of the electron, and e = — eo, eo > is its electric charge. The 
conserved total angular momentum has only a single component, namely, J z = L z + S z , 
where L z = —id/ dip and S z = <t 3 /2. 

We shall look for solutions of fl2|) in the form 



tt(t,x) 
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exp(—ieEt)ip(r, tp) 



(4) 



where e = ±1, and P > is a positive quantity. We assume the ansatz 

^(r, = e ilLp 
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(5) 
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where I is an integer number. The function ^(r,<p) is an eigenfunction of the total angular 
momentum J z with eigenvalue I + 1/2. Substituting (f|) and @ in (H), and taking into 
account of the equations 



\dr r dip . 
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we obtain 
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where a = e 2 = 1/137 is the fine structure constant. 

The exact solutions and the energy eigenvalues with eE < m corresponding to stationary 
states of the Dirac equation may be found in full analogy with the case of three space 



dimensions. We shall follow Ref. [16]. Let us look for functions / and g in the form 



/ = V^+Ee- p/2 p'- l {Q 1 + Q 2 ) , 
g = V^Ee- p / 2 p^ 1 {Q 1 - Q 2 ) , 



(8) 



where 

p = 2Ar, A = Vm 2 - E 2 , 7 = 1/2 + yj(l + 1/2) 2 - {Za) 2 . (9) 

The value of 7 is to be found by studying the behavior of the wave function at small r. From 
(0) and (D, together with the eqaulity (7 - 1/2) 2 - (ZaE/X) 2 = (I + 1/2) 2 - (Zam/X) 2 , one 
can derive the differential eqautions satisfied by Q\ and It turns out that the functions 
Qi and Q2 which rendered the solutions of flj]) finite at p = are given in terms of the 
confluent hypergeometric function F(a, b; z) as 

Q l = AF( 7 - 1/2 - (ZaE/X), 2 r , p) , 

Q 2 = BF( 1 + 1/2- (ZaE/\),2 T ,p) ■ (10) 

The constants A and B are related by 

l ZaE/X 

I + 1/2 + Zam/X v ; 

The energy eigenvalues are defined by 

1 ZaE 

7 - 5 - — = ~ n r ■ (!2) 

It is easy to show that the following values of the quantum number n r are allowed: n r = 
0, 1,2,..., if Z > 0, and n r = 1, 2, 3, . . . if I < 0. Therefore, the electron energy spectrum in 
the Coulomb field (P has the form 



E — m 

It is seen that 



1/2 

(Za) 2 

1 + 1 J 

(n r + ^/(Z + 1/2)2 -(Za) 2 ) 1 



(13) 



E = myi - (2Za) 2 (14) 

for Z = n r = 0, and E'q becomes zero at Za = 1/2, whereas in three spatial dimensions E 
equals zero at Za = 1. Thus, in two space dimensions the expression for the electron ground 
state energy in the Coulomb field of a point-charge Z\e\ no longer has a physical meaning 
at a much lower value of Za = 1/2, and the corresponding solution of the Dirac equation 
oscillates near the point r — > 0. 



III. CRITICAL CHARGE 



It is known [IB, TA\ that in three spatial dimensions the expression for the electron ground 
state energy in the Coulomb field of a point-charge Z\e\ becomes purely imaginary when 
Z > 137, and that its interpretation as electron energy no longer has a physical meaning. 
To determine the electron energy spectrum in the Coulomb field with such a charge we need 
to eliminate the singularity of the Coulomb potential of a point-charge at r = by cutting 
off the Coulomb potential at small distances. This is equivalent to taking into account of 
the nucleus size. In three space dimensions the electron energy spectrum in the Coulomb 
field regulated at small distances was first considered in [15]. With increasing Z in the 
region Z > 137, the electron energy levels in such a field were found to decrease, become 
negative, and may cross the boundary of the lower energy continuum, E = —m. The value 
of Z\e\ = Z CT \e\ at which the lowest electron energy level cross the boundary of the lower 
energy continuum is called the critical charge for the electron ground state [|12|-|14|j. If Z 
continues to grow and enters the transcritical region with Z > Z cr , the lowest electron 
energy level "sinks" into the lower energy continuum, which result in a rearrangement of 
the vacuum of the QED. This rearrangement is constrained by Pauli's exclusion principle. 
If the electron ground state at Z < Z CT is vacant, two electron-positron pairs are created; 
if it is half-occupied, one pair is created; and if it is occupied, no pairs are created. The 
Coulomb potential is repulsive for the created positrons, so they go to infinity. Hence at 
Z > Z cr a quasistationary state appears in the lower energy continuum and the new vacuum 
of QED, which corresponds to the filling of all the electron states with E < — m, has the 
total electric charge 2e [[I2|-[l4j]. Indeed, all the electron states with E < —m (the Dirac sea) 
were filled at Z < Z CI , so electrons created by the strong Coulomb field with Z > Z cr cannot 
be described by means of a convenient wave function, and the notion of charged vacuum 
was introduced to describe these states [^-pl JT^ i P^H . In terms of the new vacuum, the density 
of electric charge p(r) is classical. It is a function characterizing the spatial distribution of 
the real electric charge appearing in the new (charged) vacuum, while in terms of the old 
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(uncharged) vacuum this function should be interpreted as the probability of two electrons 
(with charge 2e) being present at a given point in space. 

We would like to see how the same system behaves in two dimensions. Let us therefore 
consider the solutions and the energy eigenvalues corresponding to stationary states of the 
Dirac equation in the Coulomb field with 2Z > 137 and find the corresponding value of 
Z cr . To find Z cr it is enough to study the energy region near the boundary of the lower 
energy continuum, —m. We shall rewrite the Dirac equation, taking account of the fact that 
eE ~ —m. Introducing functions F(r) = r/(r) and G(r) = rg(r), and eliminating G(r) 
from (0), we arrive at the equation for the function F near the boundary of the lower energy 
continuum — m in the form 

im + u _ m , + Mg^ + (gg -jo ±m F{r) = o . (15) 

We note that near the boundary of the upper energy continuum for eE « m, the function 
G{r) obeys the equation (|15D with F(r) replaced by G(r). 

Solution of (|15f) , which tends to zero at r — > oo, may be expressed by means of the 
Whittaker function (see, also \ $% ) 

F(r)~W m (2\r) , (16) 

where 



(3 = eEZa/X, u = 2^{Zaf - (I + 1/2) 2 , A = Vm 2 - E 2 . (17) 

From (7), the function G(r) at eE = —m can be obtained as 

G(r) = ±(t l+l )F-r§) . (18) 

Near the boundary of the upper energy continuum, the function G(r) is given by the Whit- 
taker function in eq. fllED with eE = m, while the function F(r) can be found from the 
relation 



Using the asymptotic representation for Whittaker function at large \z\ in the form 

W p Jz) ~ e~ z '\zf , (20) 

it is seen that the bound electron state (with \E\ < m or A > 0) is localized in the plane. 

Such a behavior of the wave function of the bound electron state may be easily understood 
if we treat ( p!q) as a one-dimensional Schrodinger-type equation which describes a particle 
with "particle energy" E' = (E 2 —m 2 )/2m in the field of the effective potential (in particular, 
for I = 0) 

Ues(r) = —eEZa/mr — (Za) 2 /2mr 2 . 

We note that the effective potential is wide enough near the boundary of the lower energy 
continuum (for behavior of the effective potential in three space dimensions see, for eg. [Q), 
and that the effective potential in two space dimensions does not contain the spin electron 
term -s(s + 1) = -3/4. 

The solution of (|15|) at eE = —m can be written in terms of the MacDonald function of 
imaginary order 

F(r) = ^K iu (\fSmZar) . (21) 



The function G(r) at eE = —m is determined by (p^8|). In the following we shall determine 
the critical value Z cr for a simple model in which the potential Ao(r) is regulated at small 
distances as follows: 

Af(r) = -Ze/r, r > R; A$(r) = -Ze/R, r < R. (22) 

In the region r < R the function F(r) obeys the equation 
d 2 F dF ( ( „ Za\ 2 9 l-l 2 



-T- + UE 



dr 2 rdr \\ R J r 2 



— -m 2 + — — )F(r) = 0. (23) 



The solution of ( |23|) is 

F(r) = r{A x JUKr) + BMkt)) , (24) 



where 

^=\/(^ + ^) 2 -™ 2 , (25) 

and J n (z) and Y n (z) are the Bessel and the Neumann functions of integer order n. 

In order for the function F(r) to be finite at the point r = we need to set B\ = 0. To 
determine the energy spectrum we need to match the solutions at the point r = R: 

(^),^-o = (^t),^ + o ' (26) 
Taking into account of the fact that R is much less than 1/m, so that k ~ Za/R, we obtain, 

for the state with I = and eE = — m, the following equation that determine (at fixed R) 

the critical charge: 

W) = ^{ i ~ v ~ z W)) • (27) 



Here X = Z cr a, v = y/AX 2 - 1, z = V8mRX, and K[ v {z) = dK iu (z)/dz. Eq. fl27|) may 
be solved numerically. As we are interested only in the critical charge corresponding to the 
ground state, we can consider small values of z. In this case, the Macdonald function with 
imaginary order K iu (z) has the following expansion: 



KiJz' 



z/smh7rz/ 



sin ( v In - + arg T(l + iv) ) (28) 



z 

1 ( . 2 



sin [v In — h arglYl + iv) + tan 1 v ) + . . . 



(29) 



4^1+^ V z 

Numerical solutions of eq. (|27|) give Z cr w 84,89 at i?m = 0.02 and 0.03, respectively. For 
comparison purpose, we recall that Z CT w 170 at Rm = 0.03 for the analogical model in 



three space dimensions [12-14 



Thus, the Dirac vacuum in two space dimensions in the presence of a strong Coulomb 
field is unstable against electron-positron production at significantly smaller values of the 
critical charge than in the case of three spatial dimensions. Another difference between 
these two cases results from the fact that electrons confined to a plane behave like a spinless 
fermion. So if the ground electron state at Z < Z cr is vacant, one pair is created; if it is 
occupied, no pairs are created. 
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IV. SUMMARY 



In this paper we present the exact solutions of the 2 + 1 -dimensional Dirac equation with 
a Coulomb field, and determine the crtical charge Z cr of a regulated Coulomb source for 
which the Dirac vacuum of the system become unstable. At Z > Z CT the lowest electron 
state of discrete spectrum is the state with n r ^ 0. So if the electron ground state at Z < Z CI 
was vacant, then at Z > Z CT an electron is created, together with a hole in the lower energy 
continuum. According to Dirac this hole is to behave as a real positive charged particle far 
from the Coulomb center. Thus, phenomena that may occur at Z > Z CT are many-particle, 
and to describe them it is necessary to apply the quantum field theory. From the point of 
view of QED, the strong Coulomb field with Z > Z CT creates a positron and changes the 
vacuum in such a way that it gains the electric charge which is exactly equal to the electron 
charge e. The spatial distribution of the electric charge appearing in the vacuum looks like 
the spatial distribution of the electron charge in the level with n r = in an atom with 
Z < Z CT . However, the density of the vacuum electric charge is a function characterizing the 
spatial distribution of the real electric charge appearing in the vacuum, while in the atom 
this function gives the probability density that the electron (with charge e) may be found 
at a given point in space. 
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